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Abstract 

One of the axioms of quantum field theory is the property of unitarity of the 
evolution operator. However, if one considers the quantum electrodynamics in the 
external field in the leading order of perturbation theory, one will find that the 
evolution transformation is a non-unitary canonical transformation of creation and 
annihilation operators. This observation was one of the arguments for the hypoth- 
esis that one should choose different representations of the canonical commutation 
relations at different moments of time in the exact quantum field theory. In this 
paper the contradiction is analyzed for the case of a simple quantum mechanical 
model being an analog of the leading order of the large-N field theory. On the one 
hand, this model is renormalized with the help of the constructive field theory meth- 
ods; the Hilbert space and unitary evolution operator are constructed. On the other 
hand, the leading order of the evolution transformation in the strong external field 
is shown to be non-unitary. Thus, unitarity of evolution in the exact theory is not in 
contradiction with non-unitarity of the approximate theory. 



e-mail: shvedov@qs.phys.msu.su 



1 



1 Introduction 



One of the main postulates of the axiomatic relativistic quantum field theory in the Wight- 
man approach |], |2| (for a review see ||) is the following. The Poincare transformations 
should be unitary operators in the Hilbert state space. For example, this postulate should 
be correct for the operator of time evolution. However, the check of this axiom for realistic 
models of QFT is difficult, since the models are usually constructed in the perturbative 
approach only. 

On the other hand, one can investigate the (3+l)-dimensional spinor QED in the strong 
external classical electromagnetic field j|, ||, |], [7|. It happens that even in the leading 
order of perturbation theory the creation and annihilation operators at different moments 
of time are related with the help of non-unitary Bogoliubov canonical transformation. 
This means that for constructing QFT in the external field, it is necessary to use different 
representations of the canonical commutation relations (CCR) at different moments of 
time. 

This observation implies that one can expect that in the non-perturbative QFT one 
should also consider different representations of CCR at different moments of time ||, 
while the time translation (evolution) is not an unitary operator in the Hilbert space but 
transformation connecting different representations of CCR. This suggestion is in contra- 
diction with the Wightman axiomatic approach. However, it is in agreement with the more 
general algebraic approach developed by Haag and Kastler |TU[ (for recent reviews of the 
algebraic approach see |TI], |12|1). 

It is shown in this paper that non-unitarity of the evolution operator in the external 
field in the leading order of perturbation theory does not contradict to the unitarity axiom 
of the "exact" theory. The simple exactly solvable model is considered in this paper. The 
Hilbert state space and unitary evolution operator are constructed with the help of the 
Bogoliubov S'-matrix approach |1| [l4| . This model can be also considered in the strong 
external field in the leading order of perturbation theory: the corresponding states are 
constructed. The evolution transformation is a non-unitary canonical transformation. 

This model corresponds to the (0+l)-dimensional "field" Q(t) interacting with infinite 
number of "fields" Qk{t)- The action of the model is 

S = I dt[L Q +f / (^-^M]-g± mk Ql (1) 



k=l 



where 



dt s 



l q = E(-i) s+1 ^(Q (s) ) 2 , Q { 

s=0 

while /ifc are some coefficients. As k —>■ oo, the set of numbers Qk tends to infinity. If 

u? a 2 

Sfcli ol — 00 > the problem of divergences arises. However, if J2T=i = oo for some m, 

k k 

the divergences can be renormalized. 

Note that the model ([!]) is a quantum-mechanical analog of the large- N theory Q(p a (p a . 
Large-iV conception is very useful in QFT: one can perform resummation of Feynman 



graphs, making use of the diagram and functional-integral techniques |TBl |T7| . The 
problem of external field (the back reaction) can be also investigated in the large- N theory 
§§■ 

Consider the theory of iV fields ip a interacting with the field <3> in the (d+l)-dimensional 
space-time. The Lagrangian of the theory is 
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Analogously to pC§ (see also |2T|, ^2], |23|, |24|]), introduce the "collective fields" being the 
operators of creation and annihilation of pairs of particles 



1 N 
kp V2N ^ k p 

where 5 k a is a creation-annihilation operator of the particle with momentum k, which 
corresponds to the field ip a . 

We will consider the states of the iV-field theory which depend on the large parameter 
N as follows 

£ / dk 1 dp 1 ...dk n dk n A+ pi ...A+ p „x£ lPl ... knPn *, (2) 

with regular as iV — > oo coefficient functions x™ and such vector \l/ that does not contain 
the particles corresponding to the fields ip a . 
Note that operators of the form 



1 N 

J dkd P y= E b£% Vkp 

V J " a=l 



multiply the norm of the state (0) by the quantity 0(N 1 ^ 2 ). Therefore, they can be 
neglected as N — » oo. In this approximation 

[^k lPl ;^k 2P2 ] - ^kik^pips + <$k lP A 2Pl )- 

Consider the free Hamiltonian H = J dkaj^Y^^i ^k°^k°> where c<j k = ^k 2 + /i 2 . If we 
consider the states of the form (||) only, it coincides with the operator 

J dkdpA^ p (uj k + uj p )A^ p , 
The operator J2a=i ¥> a (x)<^ a (x) is approximately equal to 

f dk d V (A+ p -i(k+ P )x , A- „i(k+p)xN 

(2^7 kp kp j ' 

The leading order for the Hamiltonian in 1/N is analogous to eq. ([!]): 

# = jdkd P A+ p (u k + ^ p )A kp + /dx (^n 2 (x) + f (v$) 2 (x) + ^f$ 2 

V^9 r r dk dp I A+ -i(k+p)x _i_ A- p i(k+p)x 

(27r) d J ^ j ^_ (^ kp e + /i kp e 



$fx). 



(3) 



The index A; is substituted by (k, p), the sums are substituted by integrals. Eq.([3]) can be 
also obtained from the third-quantized approach |25, 26 |. 



Investigation of the models of the type (|T]) allows us to understand the difficulties of 
the quantum field theory in the external field. 

This paper is organized as follows. In section 2 we construct the Hilbert state space and 
unitary evolution operator for the model (p]) which occurs to be renormalizable. Section 
3 deals with constructing special states of the model ([[]) which correspond to the strong 
classical external field. Time evolution in the obtained quantum theory in the external 
field is shown to be a non-unitary canonical transformation. Section 4 is devoted to the 
analysis of the paradox. 



3 



2 Construction of the model: evolution as an unitary 
transformation 

In this section the quantized model (|l|) is constructed. We show that the evolution operator 
is a well-defined unitary transformation in the Hilbert state space. 

2.1 The Bogoliubov S- matrix and unitarity of evolution 
2.1.1 Conditions on the Hamiltonian and on the state 



Models of the constructive field theory are usually formulated as follows [2~SL Instead 



of the Schrodinger equation obtained by the formal quantization procedure, 

i%= [Ho + gH^t, (4) 

where H is a free Hamiltonian, Hi is an interaction containing UV divergences, one 
considers the regularized equation. The Hamiltonian Hi is substituted by the regular 
operator Hf depending on the cutoff parameter A. At finite values of A this operator does 
not contain UV-divergences. As A — > oo, the regularized expression for the Hamiltonian 
formally tends to H\. Usually, it is necessary to add the counterterms H^(g) to the 
Hamiltonian. The evolution equation reads, 

= H^f = [H + gH? + H*(g)]*t (5) 

In the S"-matrix approach one imposes the conditions on the dependence of the counterterms 
on the cutoff parameter in order to obtain the finite ^-matrix. Within the perturbation 
framework, it is possible: the well-known Bogoliubov-Parasiuk theorem on renormalizabil- 
ity of QFT if proved. Contrary to the S'-matrix approach, in order to eliminate divergences 
from the equations of motion, it is not sufficient to impose conditions on the counterterms. 
Even in the tree Feynman graphs the Stueckelberg divergences arise []30[ when one in- 



vestigates the processes at finite time intervals (like emission of the virtual photon by a 
single electron |30fl). To eliminate the Stueckelberg divergences, one should also impose 
the conditions on the dependence on A of the initial condition for eq.(|5|). The problem of 
elimination of the Stueckelberg divergences for the leading order of semiclassical expansion 



was investigated in [[H], [!2|, |3c 

In the constructive field theory |2^, |34| one usually chooses such t-independent unitary 
transformation Ta (singular as A — > oo) that the following requirement is satisfied. Suppose 
that the initial condition for eq.(|5|) depends on A as \1/q = Ta$q , where the vector <3>q has 
a strong limit as A — > oo. Then the solution to the Cauchy problem for eq.(^) should have 
an analogous form: 

*f = f . (6) 

with regular as A — > oo vector $ A ^a^oo The vector $ 4 can be viewed as a 
"renormalized" state. The operator transforming the state <3>o to the state <3> f is regular as 
A — » oo, 

U t = s- lim U{ = hm (Ta)" 1 exp[-iH A t]T A (7) 

A^oo A^oo 

It can be viewed as a renormalized evolution operator in the model (f|). 

Note that the evolution operator exp[— iH A t] may be singular as A — >• oo while the 
renormalized operator (|7|) may be regular. 

Consider the arbitrary observable O corresponding to the operator Oa in the regularized 
theory. In the representation @ it can be written as: 

T^OaTa. (8) 
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If eq.(||) possesses the limit A — > oo, one can talk about the time-independent represen- 
tation of the observable O in the "renormalized" state space. In particular, this is a way 
to construct a time-independent non-Fock representation of the canonical commutation 
relation. 



2.1.2 The Bogoliubov approach 

To construct the operator Ta, let us use the axiomatic Bogoliubov approach [T3] (see 



also |14|) based on the conception of switching on the interaction. In this approach one 
considers the analog of the model ([|) with the time- dependent coefficient of interaction 
g t = g(t) instead of the case of the constant interaction. This generalization of the model 
seems to be a complication. However, if one considers the case of the smooth function g 
which is non-zero only on the finite time interval, the S-matrix will be regular as A — > oo, 
contrary to the evolution operator which can be viewed as S-matrix corresponding to the 
discontinuous function g being constant at t G (ii, £2) and zero at t ^ (ti, £2)- The S-matrix 
viewed as a functional on the smooth function g vanishing at sufficiently large \t\ is the 
main notion of the Bogoliubov axiomatic approach. One can determine the renormalized 
evolution operator in terms of S-matrix |14| . 



After regularization and renormalization the Bogoliubov S'-matrix takes the form 

/OO A * 

e iH °T( g ( T )Hf + H^g{-)])e- iH ^dr). (9) 
-00 

It transforms the initial condition for the equation 

*V = ^^Hi + H^[t,g{-)])e- iH ^\ (10) 

as t = —00 to the solution of this equation at t = +00, S^g)® -00 = $ +oc . Note the 
function g(t) is zero at |t| > T. 

The counterterms H^[t, #(•)]) depending on the function g and its derivatives at time 
moment t are chosen in order to make the S'-matrix regular. More precisely, for smooth 
functions g(t) the S-matrix should have a strong limit as A — > 00. 

In the interaction representation for finite values of A the evolution operator coincides 
with the Bogoliubov S-matrix @ if g{r) = g at r G [£1,^2] and g{r) = at r ^ [ti,^]- 
Namely, the substitution \1/ A = e~ lHot $ A transforms eq.(||) to the form (||), so that 

e -iH A (t 2 - tl ) = e-^S^/^.)]^^ 

where h^if) = 1 as t G (ti, t 2 ) Itit 2 (t) = as t ^ (^1,^2)- However, because of the 
Stueckelberg divergences |3(], Q the strong limit of the S-matrix as A — >• 00 for the case 
of a discontinuous function g, in general, does not exist. 

Consider the function £_(t) which switches on from to 1 at — 7\ < r < — T 2 , — T 2 < 0, 
is equal to 1 at — T 2 < r < and at r > and r < — Ti (see fig.l). 

Choose the unitary operator Ta to be the following: 

T A = SM-(-)]- (11) 

Consider also the function £+(t) = £-(— t) and operator v t of shifting the argument r: 
v t g{j) = g{r + t) . The operator U\ entering to eq.(^) takes the form: 

U^ 1 = S^[^_(.)]e-^ ot ^A[94 t2 (-)]^ otl 5A[^-(-)]- 



5 



The property S&[vtg] = e lHot Sx[g\e lHot and unitarity of the operator 5a[(?£+(-)] imply that 

U*-* = e~^St[gv t2 U-)}Sl[gv t ,A + (-)] (12) 
x5 A [^ 2 e + (-)]5 A [^ 1 i 2 (-)]^A[^ ii e-(-)]e^ 1 . 1 } 

Denote as £ tlt2 the smooth function of the form 
The operator flT2| ) can be presenter as 

Thus, the operator fll2[) is expressed via the values of the Bogoliubov S'-matrix functional 
on the smooth functions g which vanish at \t\ > T for some T. 

In the next subsection we show that the operators S\[g(-)} and S^[g(-)] have strong 
limits as A — > oo. 

This will imply that the renormalized evolution operators fl7|) exist. 

Note that two operators Ta corresponding to different functions of switching the inter- 
action ^ and ^i 2 '' lead to equivalent representations of the observables since the unitary 
operator S£[g^]S A [g^} = S+[g(^ + g^SMZ^ + &)] h as a (strong) limit as 
A -> oo. 



2.2 Construction of the Bogoliubov S'-matrix 
2.2.1 Regularization and canonical quantization 

Consider the canonical quantization of the model ([!]). Since the Lagrangian contains higher 
derivatives, the classical Hamiltonian of the model depends |35j on the coordinates Vq = Q, 
V\ = Q, V;_i = Q 1 and canonically conjugated momenta P , Pi, Pi-i, as well as 
on the coordinates and momenta q k and p k . The classical Hamiltonian function has the 
form: 

oo 

H = H q + H Q + gJ2^kq k Q, (13) 

k=l 

where 

Pk I il k q k 



tt _ y^oo [ i_k_ j 

n <l — 2^k=l I 2 + 2 

Under the canonical quantization procedure, the coordinates and momenta are associated 
with the operators V , V/_i, Pq, ...,Pi-i, q k , pk obeying the canonical commutation rela- 
tions (CCR): 

[V m ,P k ] = i5 mk , [qk,Pi]=i$ki- (14) 

Other commutators vanish. 

In order to avoid the divergences at the intermediate stages of the analysis of the model, 
introduce the regularization. The quantities are substituted by 

fi k = fi k ', k < A, = 0, k > A. 

where A is a cutoff parameter. We will show that the Bogoliubov S'-matrix is regular as 
A — > oo if the counterterm 
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is added to the Hamiltonian fll3|) . 

In the cutoffed theory one can use the Fock representation of CCR ([TJ|). The operators 
pk and <jk are presented via the creation and annihilation operators a k in the Fock space: 

They obey the following relations: 

Remind that the Fock space T is a space of sets 

* = (* ,*i(fci),*2(A;i,fc2),--.) 

of symmetric with respect to ki,...,k n functions ^ n (fei, k n ), ki,...,k n = 1,2,3,..., such 
that the series 

oo oo 

J2 E l^nih,...,^ 2 <00 (16) 
n=0 k\...kn=l 

converges. The operators a k act in the Fock space as 

n i 

(a^) n (fci, kn) = E -r^j*n-l(*li ■■; k j+U ■■; K), 
3=1 V U 

(a^*) n _i(A;i, fc n _i) = \/n^ n (k, k u fc„_i). 

The operators Vi act in the space of functions ip(V , Vj_i) from L 2 (R') as operators 
of multiplication by Vi, while p = — ^^r- 

Choose, as usual, the space T ® L 2 (R') as a state space of the composed system; the 
operators p, V and are extended as p <g> 1, V* ® 1 and 1 ® correspondingly. 

Consider eq.(|T0|) for this model. Instead of the operator e lHot (H — H )e~ lHot entering 
to the right-hand side of eq . fllQl) , where Hq — Hq + H q , consider the operator of the form 

e**{H - H 9 )e~*** = H Q + + g £ t^ 6 . (17) 

fc=l V^'fc 

Denote by C/A the evolution operator for the Hamiltonian fll7D transforming the initial 
condition at for the equation 

at t — t\ to the solution of this equation at T = t 2 , $^ = U£ t The evolution operator 
U^ t2 for eq.([T0|) is related with U^ t as 

u tit 2 — e u t 1 t 2 fi 

Let g = g(t) = gt be & function vanishing at |t| > T. Then the Bogoliubov S-matrix 
coincides with the operator U^ TT and can be expressed then via the evolution operator 
U_ TT for the Hamiltonian (0). We will show that this operator has a strong limit as 
A — > oo. This will imply that the Bogoliubov ^-matrix and the renormalized evolution 
operator for the g = const-case are well-defined as A = oo. 
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2.2.2 Renormalization of divergences in classical equations 

It is well-known || l|| that quantum theories with quadratic Hamiltonians are specified 
by their classical analogs. Consider the divergences in the classical version of the model 
(|17D . Equations of motion for the quantum Heisenberg operators coincide with the classical 
equations and have the form: 

ia~ k = 9^==Q^ kt i = S^Qe-^ (18) 



V j = V J+1 ,j = 0,l-l-V l . 1 = (-V 



a + e in k t +a - e -in k t 



-Pi = (-iy( Zj + 5z 3 )V 3 + Pj-i,j_= 1, 1 - 1; (19) 
-P = {z + 5z )Q + g J2T=i v m 

Eqs.([l9|) can be presented as 



,+„iU k t _i_ _— p -iVL k t 



E(z.qW)« + E(^ (s) ) (s) + 9 E A* a * e = o, (20) 

s=0 s=0 fc=l V z "fc 

while integration of eqs.(^) gives us the following relations: 

4(t) = of (-oo) ± i-^= r drg^^. (21) 



After substitution of formula ([H]) and integration by parts 2/ times eq. ( ^0|) is transformed 
to the following form: 



EU(^) W + £g = i^ afc( " oo)e "' t ^ ( " ooK 



/2H7 



(22) 



+ (-l) m St E^=i ft J-oo dT(g T Q T )W sm ^t^ = -(B 1 + B 2 )Q(t) 
where 

l— 1 .2 ^/2s i— 1 Js 

Check that under some choice of the counterterms B 2 + B± = 0. Note that the operator 
I?i is Hermitian and polynomial in d/dt. The degree of the polynomial is 21 — 2, so that 
B l = Y%£b j {t)4 f . Choose 

Jl-l Jl-l J21-2 J21-3 

One has 

Bt = bn-2^ + ((2l - 2)b 2l _ 2 - 6a-s)^ + ••• 

and 6 2 ;_ 3 = (/ — l)6 2 ;-2, so that the operator B\ — B^ contains the derivative d/dt in 
degrees no higher than 2l — 4. Applying to the operator B\ — B 2 an analogous procedure 
I — 2 times, one constructs such operator B 2 = B 2 + ••• + B 2 ~ 2 \ that B 2 = —B\. 
Eq.(p2|) does not contain singularities if 

Et0Ti<™- (23) 
k "fc 

Existence, uniqueness and smoothness of the solution of the Cauchy problem for this 
equation are corollaries of the general theory of integral equations (see, for example, |36| ) . 
Thus, the classical theory does not contain singularities as A — > oo, provided that the 
considered counterterms are added and condition (|23[) is satisfied. 
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2.2.3 Regularity of the Bogoliubov S-matrix 

Let us check that under condition ( p3|) the evolution operator U_ TT for the time interval 
(— T, T) is regular as A — > oo in the theory with Hamiltonian (|T7|). For simplicity of 
the notations, the index A will be omitted. Represent the operators Q = V , V s , P s via 
creation and annihilation operators B k = Vk ^[ k ; denote Bf +k = a^. The Hamiltonian 
(O) is quadratic in creation-annihilation operators: 



ij=i 



+ Re- 



consider the linear canonical Bogoliubov transformation, transforming the initial condition 
B^(—T) for the set of equations (|I8|)-(|T9|) to the solution of this set B^(T) at t — T. The 
Wick symbol of the evolution operator f/ A TT =: U\(B + , B~) : for this theory is presented 
as§]: 

u y B ' B > = 73m; (24) 

x expE£?=i \-~Bi(G2 FjfiijBj + Bi{G A l - l)ijB* + ±5*(F A G A %.B* 

where 

0j?r(T) 0j?+(T) 
^AJo" dB +(-T) ,{ A)ij dBf{-T)' 

According to the appendix, the conditions 

12ij \ (G\)ij — Gij\ 2 —>A^oo 0, YUj \ {^A-)ij ~ Fij\ 2 —>A^oo 0, 

detG\ ^a^oo detG 



(25) 



imply that the operator U^ TT has a strong limit as A — ► oo. To check condition (p5"D , it is 
necessary to show that 

(a) the I 2 - vectors of the form f^jf^ an d dQ^-T) nave strong limits as s = 0, 21 — 1. 

(b) the operators with matrices are presented as 

da H T ) (aWaMaW), da t( T ) _ n , 4(4),(5).(6K 

l^A ^A Jfcni, q i / — ^ ^A ^A Jfcm) 



9a-(-T) — "a»a^, 9<{ _ T) 
for some operators that converge as A — > oo in the norm 



\\A\\ 2 = VTrA+A. (26) 

Namely, under conditions (a) and (b) the operators Fa and G\ are evidently converge 
as A — > oo, while the determinant detG\ converges because of lemma 2 of the appendix. 
The matrices are expressed via the fundamental solution of the equation 

21 t oo 2 

5>,Ql s) ) (s) + {-l) l+1 g t / dT(g T Q T )W £ sin[fi fc (t - r)] = j t . (27) 



S = 



The solution of eq. ( B7| ) is expressed via the linear combination of the initial conditions and 
the right-hand side: 

21-1 t 

Qt= J £ Cs(t)Q (s) (-T) + / drG tT j T . 

Let the function G tT be equal to zero at t < r. Then the integration can be supposed to 
be taken from — oo to +oo. 
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It follows from eqs.(H) and (H) that 

9Q<- S HT) _ _ roo i 8 s d l (ft it _ yIm^L 
da+(-T) J ~°° <9T S dr' K^Trgr) > 



J^tiZL — A, 1V ff/r.f/T Mhe^tn fe.e'"^ 8' 9 ; / ^ 

da+(-T) ~ ° fcm H X J J ar l ar 2 ^nj+i/a V2Q^ 1/2 ^ ItfnWmSV 



(28) 



(29) 



To justify the properties (a) and (b), it is sufficient to prove the uniform convergence 
of the functions 

^1 qI ^1 ^1 

df^d?^ GTr9T ^ ^H ( ^ Cs(r)) ' d^M^ G ^ 9 ^ (30) 

at [— T, T]. Namely, the property (^) implies the convergence of vectors (p8|) . Con- 
struct operators A^K The operator transforms the sequence f m from I 2 to the 
function (A^ /)(r) = J2m ^^ttjt fm, let v4^ 2 ^ be an integral operator with the kernel 

'^Tg~r(9T 1 G T1T2 g T2 ), while the operator A^ transforms the function ip from L 2 [—T,T] to 

the sequence {A^ip) m = -ifdr ^^ f m . Choose A^ = A^ 2 \ A^ = A«, {A® f){r) = 

J2 m ^~^^TT7J~fm- Convergence of these operator in the (|26|)-norm is a corollary of the 
property ( f23"D and convergence of the functions (p0|). 

Convergence of the function J^-(g r c s (T)) is a corollary of the lemma 3 of the appendix. 
To investigate the property of convergence of the function C?(0(0 = |t|jG t17 2, represent 
eq.(^) in the form 

= (-l)'i(t - 1„), 

where (jQ is an integral operator 

((<t) /)(*) = /-oo /( r )^ r ! while is the operator with the kernel 

oo 2 

fc=l ll k 

Convergence of functions and JprJ^GW, m, n < /, is a corollary of lemma 3. 

As t > T > t, the functions G tT obey the equation Y.f=o z s-§^G tT = 0. Therefore 

r _ V-2Z-1 V / (t _ T) 2i-2m + S ^ 

^tt Z^s=0 dT^Tto ^m=[f +1] (2i-2m+s)! 2; /r>-| \ 

rtj-^ («-T)»- 1 ~ a " ^ 

Z^ m =0 JT u ' ^rto 2; (2J_l_2m)! 



The quantity -§^GTt is expressed via linear combinations of the values G rto at r G (T, T + 
5) from eq.flin]) for t = ti,...,t 2s € (T, T + 5). Therefore, the convergence of functions 
■§^-§^kGtt is a corollary of convergence of the quantity J^Gtt- 

Thus, the strong convergence of the Bogoliubov S"-matrix as A — > oo is checked for 
smooth functions g which vanish for t ^ (— T, T). Convergence of the operator S + is 
checked analogously. 
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2.3 Representations of different operators 

In this subsection we investigate what operators 0\ transform the vectors of the type T\^i 
to the vectors of the same type T^2- This means that the condition (||) is invariant under 
transformation 0\. In this case the operator (||]) is regular as A — ► oo. 

Consider the Heisenberg operators 0\ = a£(t) = e %m a\e~ im . Since the operators 
e~ lHt and Ta can be expressed via the evolution operator V tl t 2 for the theory with the 
Hamiltonian (|TT|), 

—iHt „—'iH a tjT rji jt „iHr>T 

e = e 9 Uot, Ja = tV-Toe Q 
for Heisenberg operators a£ (t) in the representation (§) one has 

T+a+(t)T A = e- lH ^ T U+ Tt e lH ^die- lH ^U. Tt e^ T 
Analogously, one has 

T+a~ k (t)T A = e- lH ^ T UU t akU-Tte lHQT e- inkt . 

TjtQ(t)T A = e~ iH ^ T U+ Tt QU. Tt e iH ^ T e- inkt . 

To investigate the regularity of the operators it is sufficient to investigate the regularity 
as A — > oo of the operators 

V+ Tt a±V_ Tt , V+ Tt QV„ Tt . (32) 

The Heisenberg equations of motion for the operators (|32"D coincide with the classical 
equations in the model (|TTD. Therefore, for operators (|32| ) one has 

(33) 

9Q(t) : - , dQ{t) A 2 f 1 dQ{t) 6(s) 
da~ m {-T) m + daU-T) m ) + to SQM(-T) V " 

It follows from the previous subsection that the operators T^a^T A and T A Q( S 'T\, s = 0, Z 
are regular as A — > oo. Equations of motion imply that the operators T A P S T A , s = 0, Z — 1, 
are singular as A — > oo. 

3 Non-unitarity of evolution in the external field 

It will be shown in this section that the evolution operator corresponding to the model Q23| ) 
in the external field in the leading order of perturbation theory may be nonunitary, since 
it is necessary to consider different representations of the canonical commutation relations 
at different moments of time. 

The quantum theory in the external field is usually constructed as follows ||. The field 
Q{t) is decomposed into two parts. The "classical" part -Q c {t) is of order 0(l/g). The 

remaining part X(t) is "quantum", 

Q(t) = -Q c (t)+X(t). (34) 



U+ Tt QU~Tt = £ 



m=l 
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The classical part of the field q k will be set to zero. Action (^4j) takes the following form: 
S = const -~JdtXj2 z s Q? s) + J dt[L x + £ | - ~p - £ ^ fc(7fc Q c ] + O(p). 

^ s=0 ' fc=l \ / fc=l 

The term of order 1/g vanishes if the "external field" Q c (t) obeys the classical equation of 
motion 

£^ = 0, 

which can be obtained from eq.(|23|) by the variation procedure as g — > 0. Neglect the terms 
of order 0(g). We obtain that the degrees of freedom corresponding to fields X and q k are 
independent. Thus, one can consider the problem of quantization of the fields qk in the 
external nonstationary classical field Q c (t). The Hamiltonian of this model has the form: 

H =£,(j + ^)+t^q k Qc(t). (35) 

k=l \ / k=l 



3.1 Fock representation: range of validity 



One can try to use the Fock representation of the canonical commutation relations (15). 
Let us investigate in what case it is possible. 

Under this choice of the representation, the Hamiltonian (|35|) takes the form: 

H = £ tl k aia- k + £ Vk^tQc + E (36) 

k k=l V Zll k 

If we choose the Wick ordering of creation and annihilation operators, the constant Eq 
vanishes. 

Consider the solution to the Schrodinger equation 

if = «, (37) 
which has the form of the coherent state 

00 

*(t) = c(t)expEa fc (t)a+]|0>, (38) 
fc=i 

being expressible via the vacuum vector |0 > of the form (1, 0, 0, ...) and complex functions 
c(t) and a k (t). Substitution of the vector ( |38"D to eq.([J7|) leads us to the relations 



ic 



00 



J2 ^gfc q + Eq 

k= i v2T^fc 



c; id k = Vt k a k + ^L- Q- 

V ^S'fc 



The divergences appearing in the multiplier c can be eliminated by the proper choice of 
the " counterterm" Eq. Investigate now the functions a k (t): 

a k (t)=a k (0)e-^ kt + Pk (t), (39) 

where 

It follows from || that expression (p8|) defines the Fock vector if 

J2\Mt)\ 2 <oo. (41) 

k 
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Integrating eq.(f40|) by parts, we obtain 

Pk(t) = ti(t)-ti(Q)e-* i >* + >i(t), (42) 

Wnere , , o 

7iW = -^/o t ^Q(^)(r)4^- 

if 



The leading in l/fi fc order is p fc (t) ~ fM=(Q{t) - Q(0)e^ fe '). Condition (P) is satisfied 



oo ,,2 



E S < oo. (43) 
fc=i 1 fc 



Thus, the evolution transformation in the external field can be viewed as an unitary oper- 
ator if the condition (031) is satisfied. 



Condition (|43[) can be also obtained as follows. Heisenberg equations of motion for the 
operators ir(a k (t)) = e lHt Tr(aj:)e~ tHt are written as 

Tij t <4{t)) = n k n(4(t)) + ^=Q(t). (44) 

Heisenberg creation and annihilation operators at different time moments are related as 

7r(af(0) = e- <n -V(a±(0))+pf(0, (45) 

where p^ = pk, pi = p%, ir(a}:(0)) = a k . According to ||, the canonical transformation 
is unitary if and only if J2k |Pfc(^)| 2 < oo. This condition is equivalent to (35). 



3.2 Different Hilbert spaces at different time moments 

If the condition (|43|) is not satisfied, one should consider non-Fock representations of CCR 
in order to construct the quantum theory. Since the choice of the non-Fock representation 



is specified by the interaction (see, for example, p7|), which depends on time in our case, 
it is necessary to consider different representations of CCR at different time moments. 
Consider this hypothesis for the model ([36]). 

One can consider the "large" linear state space C and specify the subspaces H a C C. 
The inner product is introduced on each subspace H a . The parameter a belongs to some 
set A. The operators e*^* *** defined on C transform elements of H a to elements of H a . 
The restrictions 7r Q (a^) = a fe \u a of creation and annihilation operators on the subspace 
H a specifies the a-representation of CCR. 

The evolution operator U t is defined as a set of mappings u t : A — ► A and V t a : H a — > 
TC Uta - If the initial condition belongs to H a , the state at time moment t is defined as 
V t a ^o and belongs to H Uta . 

Choose as a space L the space of analytic functionals ^(z) = *&(zi, Z2, •■■)■ The creation 
and annihilation operators have the form: 

_ d 

a+ = z k , a k = — . (46) 

Define the subset Ho G £ as follows. Consider the expansion of the functional \I/ into a 
series: 

OO 1 

* = V — — V ^ W ■ z- z- 



n=0 V ?l! i!...^ 
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with symmetric in i\...i n coefficient functions ^j™.^ . We say that \& G Hq if 

oo 

ll*l| 2 = E E l< n Uj 2 <oo. (47) 

n=0 i\...i n 



Introduce the inner product on Hq: 



(*,*)o=E*r < n*£u- 

n=0 

Note that the space Ho is isomorphic to the Fock space ||. 

Let a = (oil, OLi, ...) is a set of complex umbers. Say that \l/ G H a if the functional 

wj>(z) = e^r.i a *^K)$( 2 + a *) 
belongs to Hq. Introduce on H a the inner product: 

Eoo 
..... , , . fe =i z k a k belongs to H a in any case. This functional 

belongs to Ho if and only if a G I 2 . 

Let G H a (o). Define the mapping u t : a(0) i— > a(i) according to (|39|). Since the 

evolution equation in the representation ( ffSD has the form 

(OO Q OO £ _| 9_ \ 

E IT" + E ^-=-^-Qc + E **(*)■ 
fc=i °^ fc=1 V^fc / 

The functional $* = u> Q ,( t )\I/ < obeys the equation: 

(°° r) °° n \ 

E n***^ + E^^^c + ^oj **(*)■ 

If $° G 7io ; one has &(£) = b l ^(ze~ lUt )) for some multiplier c*. Under the appropriate 
choice of the counterterm E , $* G 7io- Thus, ^ G H a (t)- 



3.3 An algebraic approach 

The lack of the approach of the previous subsection is that it is necessary to eliminate the 
divergences from the multiplier 6* by renormalization of Eq. If one considered the density 
matrix instead of the wave function, this difficulty does not arise. The generalization of the 
density-matrix approach to systems of infinite number of degrees of freedom is the algebraic 
approach |27j which is suitable for the case when different representations of CCR arise at 
different time moments. 



3.3.1 Density matrix for systems of infinite number of degrees of freedom 

In the ci-dimensional quantum mechanics, one can use not only the wave-function language 
but also the density-matrix conception. If the system is in a pure state with the wave 
function the Blokhintsev-Wigner density matrix is determined as: 

P(P,q) = T^ Td Jd^( q -^(q+^. (48) 
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A remarkable property of the density (|4q) is that the average values of observables A 
A(p, q) can be presented in a form analogous to the classical statistical mechanics: 



< A >= (&,A&) = J dpdqA(p,q)p(p,q), 



provided that the Weyl ordering of the coordinate and momenta operators are chosen. For 
pure states, the density matrix specifies the wave function \1/ up to a multiplier. 

For the case of infinite number of degrees of freedom, the numerous difficulties with the 
divergences arise. Nevertheless, one can consider the Fourier transformation of p: 

p(a,P) = J dpdqp(p,q)e- iap -^. 

For pure states (^8|), it can be presented as 

p(a,/3) = (^,e-^- i/3 %). (49) 



The function can be used instead of the density matrix in calculations of the average 
values of observables. 

The advantage of using the function (fS|) is the possibility of generalization to the 
infinite-dimensional case. One can specify the state of the system by the average values 

z *) =< e £*(** fi * ~ z tH) > . (50) 



Consider some examples of "densities" (50). 

1. For the vacuum state 

p(z,z*) = ($(°)eSfc( z fc a fc ~ z l h l )$(°)) = e -| 

2. For the coherent state 

$ = e E fc K a fc- Q fc a fc )$(°) ) 

where a E I 2 , one has 

p(z, z*) = e£*( a i*fc- a ***-§****). (51) 

3. Suppose that the non-Fock representation n{af ) of CCR in the space TC is chosen. 
For this case, one can also use the "density" (BUI): 



p(z, z*) = ($e^k(** w (°* )-***(«*))$), $ g H. (52) 

4. As an example, consider the following "a- representation" of CCR: 

H = F, ix a (al) = al + a* k) vr a (%) = % + a k . (53) 

For the vacuum state $ = the "density" p has the form (0), but the case a G I 2 can 
be involved. 

Definition. We say that the function p(z, z*) is an a-density if it is written in the form 
for the representation (EBI). For a = 0, a-densities will be called as Fock densities. 



Statement 1. p(z,z*) is an a-density if and only if p(z, z*)eS^i'" a * 2k+n ' :2 '=' is a Fock 
density. 

Proof. The fact that p(z, z*) is an a-density means that 

p(z,z*) = ($e^'= (2fc(a fe +Q fe ) " 2: fc (a fe +afc)) $), 
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for some vector $ G T. This is equivalent to the statement that p(z, z*)eS^i(" a k 2fc+afc2 t) 
is a Fock density. 

Statement 2. Let p be an a ^-density. Then p is an a^-density if and only if a^—a^ G 

I 2 . 

Proof. According to statement 1, the function 

p(z, z*) e ^r=i(-4 1) *^+4 1) ^) = f( z , z *) = ($, e S fc fe^-^ a fe )$) (54) 
is a Fock density. Let a^ 2 -* — a*- 1 ** G Z 2 and 

One has 

Therefore, p is an a^-density. 

Let p be an a ^-density, z G Z 2 and 

Zfc 1 ^ = z k ,k < n, zj^ = 0, & > n. 

It follows from eqs. (|54|) and (^) that 

(<j> 1>e ^(«» ^ ^$0 = Efc((af - a «)4»)*-(4^-4^)4«)) 

k ( z k a k~ z k a k )$) 

It follows from the corollary 3 of lemma 1 from the Appendix that 

e £ E fc (4" )a fc-4" > s)$ 1 ) ($!, e E E fc (^ a fe- z feS 

($, e e S fc (4" )a fe-4 n) * a fc)$) ($ ) e e Sfe (2fca fe _ ^ a fe ) $) 

Corollary 3 also implies that for sufficiently small e 

e^Ek^fc-**^)^) ^ 0, 
($ )e £ E fc (^ a fc- z fc a ft)$) ^ o. 

This implies that 

E, (n) - -I- (n) * * — \ 
( fe (4 4-4 «»*)$!) ^ o, 

($ ; e £ E fe ( 2 fc a * 4 a J$) ^ 
at n > ni. Therefore, the eft-hand side of eq.(BH) tends to 



($ ; e e E fe ( 2: fe a fc- z fc a fc)<|)) 

This means that the limit 



lim Y(( 



(2) (l)s (n)* / (2)* (l)*s 



n— >oo - 1 — ' ' ' " 



(56) 



exists for all 2 G Z 2 . Choose z = i(a^ 2 ^ — a^). We see that (a^ — G I 2 . Statement is 
proved. 

We see that the notion of "density" ( |50|) is useful in order to specify states corre- 
sponding to different representations of CCR. One can even investigate the case when the 
representation is time-dependent. 
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3.3.2 Evolution of density matrix 

Let us write down the evolution equation for the average ( f)0|) for the system 
From eq. fl37| ) one has 

irho =< [ e ^ Zk ^- z t a k ) ; H] > . (57) 



for the Fock density case. Eq. (57) can be postulates for the general case as well. It follows 
from CCR that 

f( a + , a-)e^ {Zk ^~ z ^ ) = e^^" 2 ^ f(a + + a*,a~ +a). 

This implies that 



Vt k z* k a k + Vt k z k a k + Q k z k z k + p k Q — / —— 

V & 'k 



Furthermore, it follows from CCR that 







— _ < e S fc (^ a fc- z fc a fe) >=< e Efc( 2 fc a fc- z fc a fc)f a r + -z*) >■ 
dz* 2 

Substituting these relations to the commutator, we obtain the following evolution equation: 

= E (^%^ + n***^ + /^^nr J 

Substitution 

p t (z,z*) = e En=xKW**- a *( t K)/(z, z*) 
gives us eq. (|39|) on the function a k {t) and the following equation on /: 

->? = (- ntZ '*ii +n ^ f - (58) 

Statement 3. Let po be an a(0)-density. Then p t is an a(t)-density, where a{t) is given by 
eq-©- 

Proof. It is sufficient to check that the property that / is a Fock density is invariant 
under evolution fl58p. Let 

/ t = ($ t , e Er=i( a ^-w,$ t ) 

where $ t = e~Efc Cfca fe a * $ . Note that f t obeys eq.(^). Statement is proved. 
3.3.3 Time evolution of the representation 

According to statement 2, the function a k specifying the choice of the representation is 
defined up to an element from I 2 . We say that a ~ a' if a — a' e I 2 . Denote the class of 
equivalence as [a]. 

Let the condition (|23f) be satisfied. In this case the quantity j l k (t) entering to eq.( 
is an element of I 2 , so that 

[a k (t)) = [(a k (0) - Plme-^ + P l k (t)), 
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where quantities (3 l k {t) are expressed via Q(t), Q^ l ~ 2 \t). We see that in general case 
the representation of CCR at time moment t depends not only on the value of Q^ s \t) but 
also on values QS s '(0). However, there is a special case when the representation depends 
only on the derivatives of Q at the same time moment. Such a case corresponds to the 
following choice of the initial representation: 

a k (0)-/3{(0)el 2 . 

This implies that 

a k (t) - p l k (t) e l\ 

so that 

[«*(*)] = $(*)]■ (59) 

If the condition ( p3|) is satisfied at I = 2, the formula fl59|) for the representation takes the 
form: 

[«*(*)] = --£=Q(*)- (so) 

'20? 



This choice of the representation is in agreement with papers 0, 0. In these articles the 
processes in strong nonstationary electromagnetic and gravitational fields were considered. 
It was suggested to consider the representation obtained by the diagonalization procedure 
of the Hamiltonian at each time moment. 



We see also that if the condition (j23|) is not satisfied at I — 2, the prescription ( p0|) of 
Q is not valid. 



4 Explanation of the paradox 

We have constructed in section 2 the renormalized Hilbert state space and the renormalized 
evolution operator for the model (|1]) which is unitary. This implies that it is sufficient to use 
one representation of CCR. This is in agreement with the Whightman axiomatic approach. 
However, section 3 tells us that the model (0) in the strong external field in the leading 
order in g is unusual: one should choose different representations of CCR at different time 
moments. Let us discuss this paradox. 



4.1 Extraction of the c- number from the field in the cutoffed 
theory 

The procedure of extracting the c-number component of the field Q is justified within the 
Hamiltonian framework as follows [B7], ISR . Consider the regularized theory with fixed 



A. At the fixed moment of time the state of the system is specified by sets V = (Vq, ...,Vi-i) 
and Q = (Q, Q^ 1 " 1 ') and regular as g — ► vector X of the Fock space which can be 
identified with the state in the external field Q. The solution of the regularized Schrodinger 
equation (H) depends on the small parameter g as 



% = e$ St U g [P t , Q t ]X t (61) 



where the index A is omitted, 



U g [P, Q) = exp -(£(V S V S - Q {S) P S )). 



s=0 
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Substituting vector (]61|) to eq.(|]), making use of the relations 

(U g [V, Q}) + Q {s) U g [V, Q] = Q« + iQ« 



(u 9 [v,Q])+P s u g [v,Q] = P s + ±r s , m] 



we obtain that the number St is the action on the classical trajectory satisfying eqs. (|l9|) 
as g = 0: 



d-iPi 



1-1 . 



V, = Vj+i,j = 0, / - 1; Vi- X = (-1 

-Pi = {-i)Kh + Sz i) v i + p j-^3 = M^ 1 ^ ( 63 ) 
-P = (zq + Sz )Q. 

At g — > the vector X t obeys the equation: 

oo 

iX t = [H Q + H g + VkQqk + E )X t . 
fc=i 

for a c- number Eq. The operator entering to the right-hand side of this equation is presented 
as a sum of the term corresponding to the field Q and Hamiltonian (j36| ). 

Note that the classical solution Q(t) can be expressed as a linear combination of the 
initial conditions for the system (|63D: 



Q(t) = T / (a s (t)Q s + b s (t)V s ) (64) 

Analogously, the Heisenberg operators V s (t) = e lHQt V s e~ lHQt and P s (t) = e lHQt P s e~ lHQt 
obey the analog of the system (^) for operator functions. Therefore, for Vq = Q we have 

Q(t) = Y l Mt)V s + bs(t)P 8 ). (65) 

Eqs. (|64]) and ( |65"D imply that 

^ + [^, Q]s<?W W Q] = Q(*) + 0(g), (66) 
The property ( |66"D will be used in the next subsection. 

4.2 The A oo-limit 

Consider the limit A — ► oo. According to section 2, the vector \I/ A should depend on 
A according to (|6]), while $^ should be regular as A — > oo. Without loss of generality, 
consider the fixed moment of time t; index t will be omitted. Construct the state in the 
external field obeying the condition @: 

T A e^ S U g [V,Q}Y A . (67) 

The vector Y A is regular as A — > oo, Y A — >a^oo Y- Show that for finite values of A, the 
vector (|6?D is of the type ( 51 ) and therefore corresponds to the external field. Comparing 



eqs. (|67|) and (jBTj) , one obtains: 

A A = T A (V, Q)Y A , 

where 



T A (P, Q) = U+[P, Q}T A U g [V, Q}. (68) 
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The vector Y can be viewed as a renormalized state in the external field. Let us check that 
the operator T A (V, Q) is regular as A = const, g — > 0. 

It follows from eqs. (|10D and (|TTD that it transforms the initial condition for the equation 



= U ° [V ' ^^-{mi + H*[t,g(-)])e- iHot U 9 [P, Q]& A (69) 
at t — — oo to the solution of this equation as t = 0, 

$a = T A (V, Q)$ A °°. 
Making use of eq. ([T7|) , take eq.(|69|) to the form 



A = U + [Pj Q ]{gUt)Ek ^Q {t) ^ '■'+"" " 



^ s (gW(t)) 2 )f/ 9 [p,Q]$i 



(70) 



In the leading order in g one has 



,7<I)* f.+ p in k t _i_ p -in k t l-x (_-\\s 

The evolution operator V^-oo for this equation which transforms the initial state at t = — oo 
to the solution has the form 

Vt t -oo = c t exp[Y,(ak{t)at - a* k {t)a k )] (71) 

k 

where 

a${t) = -if dTUr)Q(T)-^=e^ (72) 
c* = exp( f dr[-- 2 £ t^-^ s (Q( s )(r)) 2 + ^K(r)d fc (r) - ^(r)a fe (r))]) 

J-oo 9 s= o A A 

In the leading order in g T A (V, Q) = Vo,_oo- Note that at A — > oo the operator flTIf ) 
may be singular, since the series J2k I h m A^oo & k \ 2 may diverge. 

4.3 Representations of CCR in the external field 

The average fl50|) for the semiclassical state (|67p has the form: 

< e Efc( 2 fc a fc" 2 fc a fc) >= (Y, e Sfc(^ 7r P.2( a fe)-^ 7r P,Q(ifc ))y)_ (73) 

where 

7r P , s (a±) = (T A (P, Q))+a±(T A (P, Q)). (74) 
It follows from eq. (|71|) that 

Kv,Q.{Ok ) = at + a* k , n VjQ (a k ) = a k + a k 

for ak of the form ([72]) as A = oo. We see that the function (^) is an a-density. Integrating 
eq. (|72|) by parts, we obtain that 

K] = [-^=E^Q (s) (o)^], 
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provided that the condition fl23|) is satisfied. 

This representation coincides with eq.(^). We see from the direct analysis of the exact 
quantum model that one should choose only such initial condition for the representation 



that relation ( |59| ) is satisfied 



It is interesting to note that all the representations 7Pp g are equivalent. Namely, 

W4) = U+[V, Q]n (4)U g [P, Q] (75) 

However, the operator U g is singular in g and does not possess a limit g — > 0. This means 
that representations of CCR appears to be equivalent in the exact theory and not equivalent 
in the approximate theory. 

This fact can be understood as follows. According to subsection 2.3, the structure of 
TTo(Ofc) is 

oo l—l 

7T (afc) = CLl + £ (.AkmKn + ^fcm^m) + £(CfcsVs + D ks P s ) 
m=l s=0 

for some coefficients Ak m , Bkm, Cks, Dk s - It follows from section 2 that at fixed k the 
vectors with components Ak m , Bkm, Cks, Dk s are regular as A — ► oo. Consider the limit 
A — > oo. Eqs.flTBp imply that the coefficients Ckm and Dk m are of order g, while Aki and Bki 
are of the order g 2 . It follows from eq.(|62|) that the representations (|75|) have the following 
structure: 

KV,Q( a k) = °'k + P>' + (9) ( 76 ) 
where ft = g~ l Ei=o(^ s Q (s) + D ks V s ) + 0(g) Eqs.(0) and ([7j imply that ft is 

ft = a fc (0). 

Eq. ( |7ED is in agreement with obtained in section 3. 

Thus, the representations ([75|) viewed in the leading order in g are nonequivalent if 
Sfe |ft| 2 = °° (i- e - if the condition (H3|) is not satisfied). However, in the exact theory they 
are equivalent and related with the help of the unitary operator which is singular in g. 

One can expect that analogous difficulties corresponding to extracting the classical 
component of the field arise in QED. 

This work was supported by the Russian Foundation for Basic Research, project 99- 
01-01198. 



Appendix 

In this appendix the auxiliary mathematical statements are presented. 

Lemma 1. Let U n be a sequence of unitary operators in the Fock space, and the sequence 
of their Wich symbols U n (z*,w) converges to 1 as n — > oo at arbitrary z,w G I 2 . Then 
U n — > 1 in strong sense. 

Proof. 

It follows from the Banach-Shteingaus (?) theorem (see, for example, 
sufficient to check that at some dense subset of the Fock space 

Wnf-f\\ -„->oo0. (77) 

This subset is chosen as a set of all finite linear combinations of the coherent states 

K oo 

/ = *T S m)|0>. (78) 

i=l m=l 



40 ) that it is 
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It follows from the unitarity of the operator U n that for vector (f78D the property (jT 
takes the form 

K oo 

^ a* aj exp(Y: zrzf){2 - U* n (z*, Zj ) - U(z*, Zj )) 0, 

ij=l m=l 

It is satisfied since the Wick symbol converges. 

Corollary 1. Let U n be a sequence of unitary operators corresponding to linear canon- 
ical transformations 

oo 

U-'B+Un = ^G n ) mk Bt + (F:) mk Bj:), (79) 

k=l 

the operators F n and G n with matrices {F n ) m k and (G n ) m k converge as n — > oo in the 
operator norm F n -^- n ^oo ; G n ^ n ^oo 1; for some vectors of the Fock space f and g 
(/, U n g) — > (f,g). Then U n — > 1 in a strong sense. 

Proof. According to [§], the Wick symbol of the operator U n has the form ([24]) up to 
a multiplier. It follows from lemma 1 that <0 |^ n | 0> — ► 1- Therefore, <q|'^"q> — > <o\jf\o> ■> 
< 0\U n \0 >^ 1 and U n -> 1. 

Corollary 2. Let U n be a sequence of unitary operators corresponding to linear canon- 
ical transformations (f7ty), the unitary operator U corresponds to the canonical transforma- 
tion 

U^B+U = T,kLi( G mkB£ + F^ k B^), 
F n -►„_>«, F, G n -^oo G- < 0\U n \0 >^< 0\U\0 > . 

Then U n — > 1 in a strong sense. 

Corollary 3. Let U n be a sequence of unitary operators 

oo 

U n = exp(^(4 n) a^ - 4 n) *% )), 

k=l 

while 

oo 

U = exp(J2(z k a^ - z k aj~)), 

k=l 

where G I 2 , z € I 2 , \\z^ — z\\ ^ n ^oo 0. Then U n — > 1 in a strong sense. 

To prove the corollaries, it is sufficient to consider the sequence of unitary operators 
UnU" 1 and use the statement of corollary 1. 

Lemma 2. Let A n and B n be Hilbert- Schmidt operators, which converge in the norm 
to operators A and B, while the operator 1+AB is invertible. Then det(l+A n B n ) — > n ^o 
det(l + AB). 

Proof. The property (O) Tr\XYZ\ < | \X\ \ 2 \ \Y\ \ \ \Z\ \ 2 imply that for the quantity 

_ det{l+A n B n ) 
Cn det{l+AB) 

\lnc n \ = I Jo 1 drTr([l + r(l + AB)-\A n B n - AB)]~ l 
{l+AB)- l (A n B n -AB))\ < 
max Tg (o,i) 1 1 [1 + r(l + AB)~ 1 (A n B n — Ai?)] _1 ||||(l + Ai?) _1 || 
(\\A n - A\\ 2 \\B n \\ 2 + \\A\\ 2 \\B n - B\\ 2 ) ^™ 0. 

Lemma is proved. 

Consider the sequence of the Volterra integral equations 

q n {t) = a n (t) + [ drf n (t,T)q n (r), 
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for q n at [— T, T}. Denote by G n (t,T) the Green functions for these equations which are 
defined form the relation ^ 

q n {t) = / G n (t,T)a n {r)dT 

J-T 

Lemma 3. Let the sequences of functions f n and a n uniformly converge to f and a. 
Then the sequence G n (t,r) uniformly converge to the Green function of the equation 

q(t) = a(t) + J^drf(t,r)q(r), 

and the sequence q n uniformly converge to q. 

To prove the lemma, it is sufficient to use the explicit form of the Green function which 
is obtained from the iteration procedure of [ |36|| . 
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